LSQR: An Algorithm for Sparse Linear Equations and Sparse Least Squares
نویسندگان
چکیده
منابع مشابه
An MPI-CUDA Implementation and Optimization for Parallel Sparse Equations and Least Squares (LSQR)
LSQR (Sparse Equations and Least Squares) is a widely used Krylov subspace method to solve large-scale linear systems in seismic tomography. This paper presents a parallel MPI-CUDA implementation for LSQR solver. On CUDA level, our contributions include: (1) utilize CUBLAS and CUSPARSE to compute major steps in LSQR; (2) optimize memory copy between host memory and device memory; (3) develop a ...
متن کاملSome properties of LSQR for large sparse linear least squares problems
It is well-known that many Krylov solvers for linear systems, eigenvalue problems, and singular value decomposition problems have very simple and elegant formulas for residual norms. These formulas not only allow us to further understand the methods theoretically but also can be used as cheap stopping criteria without forming approximate solutions and residuals at each step before convergence t...
متن کاملLSMR: An iterative algorithm for sparse least-squares problems
An iterative method LSMR is presented for solving linear systems Ax = b and leastsquares problems min ‖Ax−b‖2, with A being sparse or a fast linear operator. LSMR is based on the Golub-Kahan bidiagonalization process. It is analytically equivalent to the MINRES method applied to the normal equation ATAx = ATb, so that the quantities ‖Ark‖ are monotonically decreasing (where rk = b−Axk is the re...
متن کاملSparse Non-linear Least Squares Optimization for Geometric Vision
Several estimation problems in vision involve the minimization of cumulative geometric error using non-linear least-squares fitting. Typically, this error is characterized by the lack of interdependence among certain subgroups of the parameters to be estimated, which leads to minimization problems possessing a sparse structure. Taking advantage of this sparseness during minimization is known to...
متن کاملA matrix LSQR algorithm for solving constrained linear operator equations
In this work, an iterative method based on a matrix form of LSQR algorithm is constructed for solving the linear operator equation $mathcal{A}(X)=B$ and the minimum Frobenius norm residual problem $||mathcal{A}(X)-B||_F$ where $Xin mathcal{S}:={Xin textsf{R}^{ntimes n}~|~X=mathcal{G}(X)}$, $mathcal{F}$ is the linear operator from $textsf{R}^{ntimes n}$ onto $textsf{R}^{rtimes s}$, $ma...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: ACM Transactions on Mathematical Software
سال: 1982
ISSN: 0098-3500,1557-7295
DOI: 10.1145/355984.355989